Cooperative enhancement of channeling of emission from atoms into a nanofiber 
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We show the possibility of directional guided superradiance from a string of atoms separated by 
one or several wavelengths in a line parallel to the axis of a nanofiber. We find that the rate and 
efficiency of channeling of emission from the atoms into the fiber are cooperatively enhanced by the 
guided modes. 
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Coupling of light to subwavelength structures and its 
control pose one of the greatest challenges of recent re- 
search []], 0, 0, 0, IE 0] • Strong coupling in a supercon- 
ducting circuit at microwave frequencies has been ob- 
served [2j. Chang et al. have proposed a technique 
that enables strong coherent coupling between individual 
emitters and guided plasmon excitations in conducting 
nanostructures In the case of dielectric waveguides, 
it has been shown that a significant fraction of emission 
from a single atom can be channeled into a nanofiber 
0) H, @|- The formation of single-atom trapping sites 
on a nanofiber surface without any external field and 
the coupling of single photons to a nanofiber have been 
demonstrated . The cooperation of two distant atoms 
via a nanofiber has been discussed 0] ■ It has been shown 
that, at large distances between the atoms, a substantial 
energy exchange can survive due to the guided modes 
In this paper, we show the possibility of a directional 
guided superradiant emission process that can enhance 
the rate and efficiency of channeling of emission from a 
string of distant atoms into a nanofiber. 

Before we proceed, we note that superradiance is a 
problem of fundamental interest [1| . Despite a great deal 
of research 0, [l(| , certain aspects of the problem are still 
not well understood. An example is the mode selection 
in the directional emission from an extended sample of 
atoms. The main difficulty is due to the fact that the 
collective process involves a huge number of degrees of 
freedom, associated with many atoms and a continuum 
of field modes. Recently, the angular distribution of emis- 
sion from a spatially extended array of atoms in free space 
has been treated by the quantum trajectory approach 
pl| . The dynamic mode selection has been studied [l2T |. 
Superradiant conversion of atomic spin gratings into sin- 
gle photons in an optical cavity has been demonstrated 

Consider N identical two-level atoms interacting with 
the quantum electromagnetic field in the vicinity of a 
nanofiber (see Fig. [1}. The fiber has a cylindrical sil- 
ica core of radius a and refractive index n\ = 1.45 and 
an infinite vacuum clad of refractive index ria = 1. Wc 
assume that the atomic frequency luq is well below the 
cutoff frequency of the fiber, so the single- mode condi- 
tion is satisfied for this frequency. In view of the very 
low losses of silica in the wavelength range of interest, we 
neglect material absorption. The atoms are located at 



points (rj, tpj, Zj), where j — 1,2, . . . ,N labels the atoms 
and (r, (p, z) is the cylindrical coordinates with z being 
the axis of the fiber. We assume that the field is initially 
in the vacuum state. The field can be decomposed into 
the contributions from the guided and radiation modes, 
whose quantum expressions are given in Ref. Q ■ 
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FIG. 1: String of atoms in the vicinity of a nanofiber. 

Assume that the characteristic atomic lifetime is large 
as compared to the optical period 2n/ui and to the 
light propagation time between two different atoms. The 
master equation for the reduced density operator p of 
the atomic system in the electric-dipole, rotating-wave, 
and Born-Markov approximations has been previously 
derived @, S Gj]> Ell • I n the interaction picture, it reads 



1 N 



j ij (2cr j pal - ajo-jp - paja-j). 
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Here, aj and en are the pseudospin operators that 
describe the downward and upward transitions of the 
atoms. The coefficients 7^ = 7^ gmdod ^ + 7^ ad \ with 
i, j = 1, 2, . . . , N, characterize the collective spontaneous 
emission process, where 7j? md and 7y d ' arc the con- 
tributions from the guided and radiation modes, respec- 
tively [7[ ■ In Eq. ([T]) , we have neglected the propagation 
effects and the dipole-dipole interactions. For these ap- 
proximations to be valid, the sample must not be too 
large and the distances between the atoms must not be 
too small. We note that Eq. has the same form as 
that for atoms in free space [ll|, [lj] . 

We introduce the total emission intensity / = 
^2huj a (h a }, the intensity of emission into the guided 
modes /guided = S guided huj a (n a ), and the intensity of 
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emission into the radiation modes / r ad = Erad husaiha). 
Here, uj a and (n a ) are the frequency and mean num- 
ber, respectively, of photons in mode a. We find 

/ = fnOoJ2i J lij(< J i (J j) = ^guided +^rad, where /guided = 

^oE 4 j7l| Uldcd) (4 cr j} and / ra d = hua Ey lij 
We note that / = —TiujqP, where P = Ej( cr j cr j) 1S trie 
total population of the excited levels of the atoms. The 
total energy emitted from the atoms is U = J °° I(t) dt = 

^guided + C^rad, where t/guided = Jo°° guided (*) dt and 
f/rad = f Imd{t)dt are the energies emitted into the 
guided and radiation modes, respectively. The fractions 
of energy emitted into the guided and radiation modes 
are given by /guided = U m ided/U and / rad = U Ta d/U = 
1 - /guided, respectively. 

The diagonal coefficients jjj describe the spontaneous 
decay of individual atoms. The off-diagonal coefficients 
7jj', with the convention j f, characterize the energy 
transfer between two atoms. According to Ref. Q, the 

contribution 7 J y Ulded ' ) of the guided modes to the trans- 
fer rate is periodic in the z direction with the period 
Xf = 27r//3o, where @q is the longitudinal propagation 
constant of the guided modes at the atomic frequency 



(rad) 



of the radiation 



ujq. Meanwhile, the contribution 7 - ■, 
modes reduces to zero with increasing interatomic dis- 
tance \zj — Zji |. Therefore, in the limit of large \zj — zy , 
the transfer coefficient 7^' is mainly determined by the 

contribution 7j|, ulded ' of the guided modes and is almost 
periodic with the spatial period Xp. 

We now assume that the atoms are aligned along a 
line parallel to the fiber axis, with relatively large atomic 
separations being equal to integer multiples of the longi- 
tudinal wavelength A^. In other words, we assume that 



ro, <fj = const = ipo, and Zj+± 



1jX F , 



Tj = const 

with qj being nonzero, positive integer numbers. In ad- 
dition, we assume that the dipoles of the atoms are ori- 
ented in the same direction. Under these conditions, the 
guided transfer coefficients 7 J y Uldcd ^ achieve their axial 
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while, due to the large separations between the atoms, 
the radiative (unguided) transfer coefficients 7^-/ are 
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small. In this case, we have jjj — 7^ 

j r^> (guided) 

7 = 7guided + 7rad and Tjj/ = 7V?, 

7j?/ ded) = 7guided- Here, 7 gU ided and 7 rad are the rates 
of decay into the guided and radiation modes, respec- 
tively. They do not depend on the axial coordinate z 
of the atoms, but increase with decreasing atom-surface 
distance r — a. We display in Fig. [3] the rates 7 gu ided, 
7rad, and 7 = 7 gu ided + 7rad, calculated for a radially 
oriented dipole with the cesium D2-lme transition wave- 
length Aq = 852 nm in the presence of a fiber with radius 
a = 200 nm. In particular, for the atom-surface dis- 
tance r — a = 100 nm, we obtain 7 gu ided = O.2670 and 
7rad = I.O67Q. Here, 70 is the atomic natural linewidth, 
whose magnitude is about 5.3 MHz in the case of the 
cesium D2 line. We note that the creation of regular 



strings of atoms in a standing wave optical dipole trap 
has been demonstrated The formation of single- 

atom trapping sites on a nanofiber surface without any 
external field has been reported Q. Superradiance of 
lines of atoms in free space has also been studied [ill ]. 
We emphasize that the transfer coefficients 7^/ for dis- 
tant atoms considered here are mainly due to the guided 
modes and are substant ially larger than those for distant 
atoms in free space [ll|, Hif . 

For the string of atoms described above, we find that 
the intensity of emission into the radiation modes is 
^rad = ^wo7rad-P and, hence, the intensity of emission 
into the guided modes is 



dP 

^guided = -1U0Q ( — + 7rad-P 



(2) 



Meanwhile, the total excited-state population P is gov- 
erned by the equation 



dP 
~dt 



-7-P-7guided5m Cr j')- 



(3) 
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FIG. 2: Total decay rate 7 (solid line) and the contributions 
7guided (dashed line) and 7 ra d (dotted line) from the guided 
and radiation modes, respectively, as functions of the atom- 
surface distance r — a. The fiber radius a = 200 nm and the 
atomic wavelength A u = 852 nm are used. The dipole of the 
atom is radially oriented. 

We examine two cases of initial atomic states. First 
we consider the case where the atomic system is ini- 
tially prepared in the symmetric one-excitation state 
|1) = iV-V2 Si |i.). Here, |1 3 ) = |e 3 -) ® \g f )_ is 

the product state in which only atom j is excited, with 
\ej) and \gj) being the excited and ground states, respec- 
tively, of atom j. The state |1) is the first excited Dickc 
state [§[ and is an entangled state, which is of great in- 
terest in quantum information and quantum computation 
0,113 ■ We introduce the notation |0) = J7, \9j) f° r tne 
state in which all the atoms are in their ground states. 
We find that the N two-level atoms prepared in the state 
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|1) act like a single effective two-level system, with the 
upper level |1) and the lower level |0). We obtain from 
Eq. (JlJ) the solution pn = e~ rt , poo = 1 — e~ rt , and 
Pw = Poi = 0, with the collective decay rate 



r = 7r 



ad 



^7g 



ided ■ 



(4) 



This rate is enhanced [8( by the cooperativity of the 
atoms via the guided modes. The above solution yields 
the total excited-state population P = e~ r * and the in- 
tensity of emission into the guided modes 



^guided — ?iWo-/V7 gu idede 



(5) 



Hence the energy emitted into the guided modes is 
^guided = SwoA f 7 gu ided/r- Meanwhile, the total emit- 
ted energy is U = Tujq. Consequently, the fraction of 
energy emitted into the guided modes is 



/guide 



-^7guided 



7rad + A r 7g U ided 



(6) 



It is clear that / gu ided increases with increasing atom 
number N and that / gu ided — * 1 in the limit N — » oo. 
Thus the efficiency of channeling of emission from the 
atoms into the fiber is cooperatively enhanced. We use 
Eq. © to calculate / gu ided as a function of N for the 
parameters of Fig. [2J and display the results in Fig. 
ill For N = 100 and r - a < 200 nm, we obtain 
/guided > 0.92 (see the endpoints of the curves). In par- 
ticular, for N = 100 and r — a = 100 nm, the factor 
/guided reaches the value 0.96 (see the endpoint of the 
dashed curve). Such a high efficiency indicates that the 
single photon emitted from the atoms is almost entirely 
directed into the guided modes. A very similar result 
has been obtained for the superradiance of atoms in an 
optical cavity [l3j]. Indeed, in terms of the single-atom 
cooperativity parameter rj = 7 gU ided/7rad) the channeling 
efficiency / gu ided given by Eq. © coincides with the suc- 
cess probability V = Nr)/(l + Nrj) for conversion in the 
cavity case [13| ■ Such a coincidence is due to the fact that 
the nanofiber mode and the cavity mode have many com- 
mon features. We note that, at the distance of 100 nm 
from the surface of the 200-nm-radius fiber, the single- 
atom cooperativity parameter is r\ — 0.25. This value is 
substantially larger than the value r\ — 6.9 x 10~ 3 for a 
moderate- finesse cavity (T3 |. 

We now consider the case where all the atoms are ini- 
tially prepared in the same coherent superposition state, 
that is, the initial state of the atoms is the product state 
= n,'( cos § \ e j) + e ^ sm f Iffj))- Such a state can be 
prepared by using a plane- wave optical pulse to excite the 
atoms. We solve Eq. (JTJ numerically for N — 10 atoms, 
use this solution to calculate the intensity of emission into 
the guided modes / gu ided, and show the results in Fig. [4] 
(solid curves). For comparison, we also show the results 
for a single atom (dashed curves). The comparison be- 
tween the solid and dashed curves shows that the energy 
emitted into the guided modes, determined by the area 
under the intensity curve, is enhanced by the collective 




FIG. 3: Fraction /guided of energy emitted from atoms ini- 
tially prepared in the symmetric one-excitation state into the 
guided modes as a function of the atom number N. The 
atom-surface distance is r — a — 200 nm (solid line), 100 nm 
(dashed line), and (dotted line). The parameters used are 
as in Fig. 



effect. Typical features of superradiance, such as the in- 
crease of the emission rate, the occurrence of a local peak, 
and the enhancement of the peak intensity, are observed. 
However, they are rather weak in the case of Fig. 01 The 
reason is that the number of atoms iV is not large and the 
transfer rate 7 gu ided is small compared to the radiative 
decay rate 7 ra d- We expect that the use of larger values 
for N would lead to more dramatic effects. However, it 
is difficult to find the exact numerical solution to Eq. (JTJ) 
for large N because the number of variables is large. 

In order to get insight into the case of large N, we make 
an approximation for the last term in Eq. (J3J). For the 
initial product state |*}, we have (<j)<Tj>) = P{N-P)/N 2 
for every pair j =^ f. We assume that this relation is 
valid for the whole emission process [l(| [III]. Such an 
assumption is reasonable under the condition N 3> N — 
Pa > 1 [9] . With this assumption, Eq. |J3J) yields 



dP 
~dt 



-P[7+(l-7V- 1 )7 guidod (7V-P)]. 



(7) 



The solution to the above equation, subject to the initial 
condition P(0) = P , is P = N(n + l)/[n + e r(t+t ^}, 
where k= {N-l)-y guidcd /-f and t a =t\o[(k+1){N/Pq)- 
k], with t = r . The intensity of emission into the 
guided modes is 
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K + e r(t+t a ) 

p T(t+t 

7(« + l)- 



T(t+t a 



7rad 



(8) 



If t a < t p , where t p = rln{(l - jV -1 )[2 + (N - 
2)(7 g uided/7)]}; then the intensity I gu ided{t) has a local 
peak with the height P^ cd = hcj j gui dcdN 3 /[4(N - 1)] 
at the time t = t max = t p — t a . Otherwise, the function 
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FIG. 4: Intensity of emission from N = 10 atoms, initially 
prepared in a product state, into the guided modes as a func- 
tion of time (solid lines). The angle 9 for the initial product 
state is (a) 9 = (full excitation) and (b) 9 — n/2 (one-half 
excitation). The distance from the atoms to the fiber surface 
is r — a — 100 nm. The intensity is normalized to NIo, where 
lo = hioo'yo- The time is normalized to to = "fg 1 ■ The pa- 
rameters used are as in Fig. [2] For comparison, we also show 
the results for the case of a single atom (dashed lines). 



traverse the sample in a time shorter than the charac- 
teristic time scale of the collective decay. Therefore, 
the length L of the atomic string in our model is lim- 
ited by the condition L -C Lq, where Lq = c/T is 
the cooperativity length. When we take N — 100 and 
7o = 5.3 MHz, and use the parameters 7 gu idod = O.2670 
and 7 ra d — I.O670, obtained in Fig. [2] for the atom- 
surface distance r — a = 100 nm, we find Lq = 33 cm. 
For L > Lq, the collective effects can still survive but the 
propagation effects must be included. 
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^guided (i) monotonically decreases from its initial value 
/guided(O) = ft Wo P 7 guidod [l + (N- 1)(1 - Po/N)}. Very 
similar features have been obtained in the case of atoms 
in free space [HI EH ■ 

It follows from Eq. © that the energy emit- 
ted into the guided modes is f7 gu ided = huoPo{l — 
N p™ d ^ In i + (i^pp/jv)« }• Meanwhile, the total emitted 
energy is U = TiluqPq. Hence the fraction of energy emit- 
ted into the guided modes is 



/guide 
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N 7rad 1 
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Po 7 « l + (l-P /iV> 
Under the condition Po = N, we have 

7rad ln(« + 1) 



/"guided 1 



7 



(9) 
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Equation (flO|) together with the expression k — (N — 
l)7guided/7 indicate that / gu idod increases with increasing 
N and that / gu ;dcd — > 1 in the limit N — > 00 (see Fig. [5]). 

Since the propagation effects are neglected in our 
model, the above results are valid only if a photon can 



FIG. 5: Fraction /guided of energy emitted from atoms ini- 
tially prepared in a product state into the guided modes as 
a function of the atom number N. The angle 9 for the 
initial product state is (a) 6 = (full excitation) and (b) 
9 = 7r/2 (one-half excitation). The atom-surface distance is 
r — a = 200 nm (solid line), 100 nm (dashed line), and 
(dotted line). The parameters used are as in Fig. [2] 



In conclusion, we have shown the possibility of direc- 
tional guided superradiance from a string of atoms that 
are separated by one or several wavelengths in a line par- 
allel to the axis of a nanofiber. The rate of emission is en- 
hanced by the cooperativity of the atoms via the guided 
modes. The efficiency of channeling of emission into the 
guided modes increases with increasing atom number and 
approaches unity in the limit of large numbers of atoms. 
In particular, for a string of 100 atoms prepared in the 
symmetric one-excitation state at the 100-nm distance 
from the surface of a 200-nm-radius fiber, the channeling 
efficiency can be as high as 96%. 

This work was carried out under the 21st Century COE 
program on "Coherent Optical Science." 
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